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ABSTRACT 

In this thesis models of MFSK digital communications 
systems were developed using Matlab Simulink and 
Communications Toolbox. The models were employed to verify 
MFSK performance in the presence of additive noise and 
predict MFSK performance for additive noise and co-channel 
interference. Results are presented for bit-error rate as 
functions of the signal-to-noise and signal-to-interf erence 
power ratios. The results for coherent detection of MFSK 
(M=2,4,8) in the presence of additive white Gaussian noise 
show excellent agreement with the theory. On the other hand, 
simulation results for the probability of bit error for non- 
coherent detection of MFSK differ (-21% average) from the 
theory suggesting a possible 'systematic' error in the 
Communications Toolbox implementation of the non-coherent 
MFSK detection. 



v 




vi 






DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
TABLE OF CONTENTS^^^^Y CA 93943-5101 

I. INTRODUCTION 1 

II. SIMULINK AND COMMUNICATIONS TOOLBOX 3 

III. MODEL OF A MFSK PASSBAND DIGITAL COMMUNICATION AND 

ITS BASEBAND EQUIVALENT 19 

A. BASEBAND EQUIVALENCE FOR PASSBAND SIGNALS AND 

SYSTEMS 19 

B. MONTE CARLO SIMULATION 22 

C. M-ARY FREQUENCY- SHIFT KEYING (MFSK) 25 

1. Bit Error Probability for Coherent MFSK 26 

2. Bit Error Probability for Non-Coherent 

MFSK 28 

D. MFSK MODELS AND BLOCK ANALYSIS SIMULATION 29 

1. Coherent MFSK Model (Noise Case) 29 

2. Block Analysis 30 

3. Coherent MFSK Model for MFSK with 

Interference and Additive Noise 39 

IV. SIMULATION, ANALYSIS , AND PERFORMANCE VERIFICATION 

FOR WIDEBAND NOISE 45 

A. PERFORMANCE VERIFICATION FOR WIDEBAND NOISE 45 

1. Results for Coherent Binary FSK (BFSK) 46 

2. Results for Coherent 4FSK 50 

3. Results for Coherent FSK. 55 

4. Results for Noncoherent Binary FSK (BFSK) ...60 

5. Results for Noncoherent 4FSK 69 

6. Results for Noncoherent 8FSK 74 

B. OBSERVATIONS 79 

V. SIMULATION AND PERFORMANCE ANALYSIS FOR CO-CHANNEL 

INTERFERENCE 83 

vii 



A. INTERFERENCE IN DIGITAL COMMUNICATION SYSTEMS 83 

B. PERFORMANCE OF MFSK DIGITAL COMMUNICATION SYSTEMS 

IN THE PRESENCE OF ADDITIVE GAUSSIAN NOISE AND CO- 
CHANNEL INTERFERENCE 84 

1. Results For Coherent Binary FSK (BFSK) 85 

2. Results For Coherent 4 FSK 8 9 

3. Results For Coherent 8FSK 92 

4. Results For Noncoherent Binary FSK (BFSK) ...96 

5. Results For Noncoherent 4FSK 99 

6. Results For Noncoherent 8FSK 102 

C. OBSERVATIONS 106 

VI. SUMMARY AND CONCLUSIONS Ill 

APPENDIX A 113 

APPENDIX B 115 

APPENDIX C 117 

APPENDIX D 119 

APPENDIX E 123 

LIST OF REFERENCES 127 

INITIAL DISTRIBUTION LIST 129 



ACKNOWLEDGMENT 



Numerous individuals assisted me in completing this 
thesis. Many thanks to Dr. Jovan Lebaric for his technical 
support, clear explanation and supervising. I would also 
like to thank Prof. Clark Robertson and Prof. David Jenn for 
their advice and help. 

Finally, I would like to thank my wife Antigoni who 
always supports me with love and encouragement, my parents 
and my country who gave me the means to acquire my 
knowledge . 



IX 



I. 



INTRODUCTION 



A. OBJECTIVES 

In many military scenarios a large number of different 
communication systems operate simultaneously from the same 
site or platform, ideally with as little mutual interference 
as possible under the constraints of available frequency 
allocations and the physical separations of the antennas. It 
is important to predict the effects of co-channel 

interference such that measures may be taken to reduce the 
interference to a tolerable level. 

The objective of the thesis is to model an M-ary 
frequency-shift keying (MFSK) communication system in the 
time domain using Simulink and MATLAB/SIMULINK 

Communications Toolbox. MFSK performance is verified for 
additive white Gaussian noise (AWGN) . In addition, MFSK 
performance for simultaneous AWGN and co-channel 
interference is determined. This is accomplished by 
obtaining estimates via computer simulation of the bit error 
rates under these conditions. 

B . METHOD 

Using Simulink and the Communications Toolbox allows 
one to concentrate on system modeling rather than on the 
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details of code development in Matlab or some other computer 
language. The block-diagram based SIMULINK time domain 
modeling is a convenient tool for the visualization of 
communication signals at various stages of transmitters and 
receivers and for "Monte Carlo" type simulations of their 
performance. The results of Monte Carlo type computer 
simulation were used to establish the bit error rates, under 
realistic conditions of noise and interference and for 
different values of transmitter/receiver and channel 
parameters such as the power ratios of the desired and 
interfering signals. 

In Chapter II the Simulink and Communications Toolbox 
is described. Chapter III discusses the models used for the 
simulation. The simulation results for coherent and 
noncoherent MFSK (M=2, M=4, M=8) in the presence of AWGN are 
verified in Chapter IV, and in Chapter V the results are 
discussed for the case involving the presence of both AWGN 
and interference. Finally, Chapter VI presents a summary and 
conclusions . 
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II. SIMULINK AND COMMUNICATIONS TOOLBOX 



A. INTRODUCTION 

Simulink is a program for modeling and simulating 
dynamic systems including linear, nonlinear, discrete-time, 
continuous-time, and hybrid systems. As an extension to 
Matlab, this environment adds many features specific to 
dynamic systems while retaining all of Matlab' s general- 
purpose functionality. Model definition and model analysis 
are the two steps in a Simulink computer simulation. In 
practice, these two steps are often performed interactively 
as the model is created and modified in order to achieve the 
desired behavior. 

Simulink represents systems as block diagrams. Using 
mouse driven commands, one can create models whose 
parameters can be edited or changed using the keyboard. 
Also, operations can be clicked-and-dragged, and results can 
be displayed in real time during a simulation. The models 
can be analyzed either by choosing options from the Simulink 
menus or by entering commands in Matlab' s command window. 
The results, which can be viewed while the simulation is 
running, are made available in the Matlab workspace after a 
simulation has been completed. Since Simulink is built on 
top of the Matlab environment, it has many capabilities. 
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which can be further enhanced with application-specific 
toolboxes. Generaly speaking, Simulink combines the power 
and ease-of-use of an application package with the 
flexibility and extensibility of a language [Ref 2,1]. 

B. TIME DOMAIN SIMULATION 

In order to evaluate the time domain simulation, the 
terms "modules" and "single-point format" need to be 
explained first. The modules of a simulation program can be 
regarded as subroutines that operate on the data, which is 
provided as input. The program calls the first module, which 
processes the data and returns the results. Then it calls 
the second module to process the data delivered from the 
first module and so on. This continues until all modules 
have been called, upon which the program repeats from the 
beginning with new input data. The term "single-point" 
refers to the amount of data that can be processed on each 
call of a module. When only one data point (only one time 
instant) can be processed or one data point is passed for 
computation from module to module, such simulation format is 
called "single-point data." On the other hand, if a whole 
block of N data points (N time instants) is processed on 
each call, the format is called "block data." This means 
that for an equal amount of total data, single-point format 
calculations are substantially more time consuming than the 
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block data operations since single point data calculations 
incur at least N times as much calling overhead as block 
data simulations for the same (N) number of time instants. 
The single-point format (which is the only alternative when 
a system includes feedback loops) processes only time-domain 
signals (periodic or aperiodic) . Usually, single-point data 
is regarded as one sample from a long, aperiodic signal, and 
thus a single-point simulation will exhibit the same 
transient behavior as the system being simulated. Also, 
single-point simulations can be arbitrarily long in duration 
without special provisions. This makes the simulation more 
convenient when collecting data for performance statistics 
[Ref 6] . 

C. CONSTRUCTING A SIMPLE MODEL 

Simulink uses block diagrams to represent dynamic 
systems. Defining a system is much like drawing a block 
diagram. The "building blocks" can be copied from various 
block libraries, which can be either built-in (Simulink or 
Communications Toolbox libraries) or user-created. The 
built-in Simulink is organized into sub-libraries according 
to block function or nature. In order to assemble a system, 
blocks can be copied from an existing library and connected 
using the mouse and new block parameter values inserted in 
place of the default values, thus saving time when building 
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new models. As an example, the process of setting up a 
simulation of a very simple system is outlined. Suppose that 
the output of a pulse generator is to be observed on an 
oscilloscope. First, the command Simulink is typed at the 
Mat lab command window prompt in order to start Simulink and 
access the built-in Simulink libraries. A double-click on 
these basic libraries opens another window containing the 
related blocks. Figure 1 shows the basic Simulink libraries 
(which show upon typing the command Simulink ) . Figure 2 
shows some of the blocks available in the Sources library. 




Figure 1. The Standard Simulink Block Library. 
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Figure 2. Blocks Available in the Sources Block Library. 

To build the simple simulation the Sources sub-library 
needs to be opened. This displays yet another Simulink 
window from which the discrete pulse generator can be copied 
to a "New File" window. Similarly, by opening the Sinks sub- 
library, the oscilloscope block can be copied and placed 
into the new system window. Figure 3 shows the discrete 
pulse generator placed into a new window. Figure 4 shows the 
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oscilloscope block placed into the new window and finally 
Figure 5 shows the finished model, with the source connected 
to the oscilloscope. 



, exaiiibte I d 






PSEfSOI 







Dlscrtt* Puls* 



Figure 3. Placing the Discrete Pulse Generator Block. 




Figure 4. Placing the Oscilloscope Block. 
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Figure 5. Completed Simulink Model. 

All copied blocks maintain the original default 
internal parameters until modified by the user. Therefore, 
by simply connecting the sine wave "source" block to the 
oscilloscope "sink" block and selecting the parameters of 
the source, the system is ready to perform a simulation [Ref 
1 , 2 ] . 



9 





D . COMMUNICATIONS TOOLBOX 

The Communications Toolbox is designed for use with 
Matlab and Simulink and represents a collection of 
computational functions and simulation blocks for research, 
development, system design analysis, and simulation in the 
communications area. The toolbox contains ready-to-use 
functions and blocks, which can be easily modified by the 
user to implement the required schemes, methods, and 
algorithms. Simulink blocks and Matlab functions accelerate 
the design process by helping users to rapidly develop and 
analyze different system designs. By using the 
Communications Toolbox (whose functions can be either called 



directly from 


the Matlab 


workspace 


and .m 


files or 


implemented 


as 


specialized 


Simulink 


blocks) time-domain 


simulation 


of 


various communication 


systems 


can be 


implemented. 


The 


toolbox adds 


a variety 


of Matlab 


functions 



and Simulink blocks that are very useful for communication 
system simulation. These include: 

• Signal Generators (data source) 

• Source Coding/Decoding 

• Error-control Coding 

• Baseband and Passband Modulation/Demodulation 

• Transmitting and Receiving Filters 
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• Baseband and Passband Channel Models 



• Multiple Access 

• Synchronization 

• Galois fields calculations. [Ref 1] 

E. CONSTRUCTING A COMMUNICATIONS SYSTEM MODEL 

A simple communication system contains these basic 
components : 

• information source, 

• the transmitter, 

• the channel, 

• the receiver, and 

• information "sink." 

Suppose that the performance of a coherent MFSK 

communication in the presence of additive white Gaussian 
noise needs to be determined. The Communications Toolbox is 
started by typing the command commlib at the Matlab command 
line prompt. This gives access to the "top level" 
communications library shown in Figure 6. The transmiting 
part of a communication system is comprised by the upper row 
of blocks while the receiving part is comprised of the 
bottom row of blocks. 
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Figure 6. The Communications Toolbox Simulink Block Library. 



Since in most cases the receiving computation is 
exactly the reverse of the transmitting computation, a 
double click on any of the blocks in this library opens a 
sub-library that includes both the transmitting and the 
corresponding receiving blocks. A sub-library typically 
contains function blocks used to build the models and demo 
blocks (cyan colored) that are used to demonstrate a block 
operation. Finally, by opening the source/sink library 
various types of sources become available. In order to build 



12 





an MFSK system the "Sampled Read from Workspace" block is 
copied from the Source sub-libarary and placed into a new 
Simulink window. The "Sampled Read from Workspace" block 
reads a row of data from the Matlab workspace at every 
sampling point. Similarly, by opening the modulation/demo- 
dulation library (Figure 7) and double clicking on the 
Digital Modulation/ Demodulation sub-library, the MFSK Mod 
Baseband and the Coherent MFSK Demod blocks can also be 
copied and placed into the new window. 



n Library: com modu BUB 



£te Edit Simulation ;• 

Modulation, ©emoduiation Library 



passband baseband 

Analog Mo/Dem: 



Analog 
mo/d em 



Analog 
mo/dern 
cmplx env 



Digital Mo/Dem: 



Digital 

mo/dem 

M 



Digital 
mo/dem 
cmplx env 



Separated Version for Digital Mo/Dem: 



Digital 




Digital 


mo/dem 




mo/dem-ce 


map/demap 




map/demap 



Note: 



Digital Demodulation = demodulation+demap 
Digital Modulation = map+ modulation 



Figure 7. The Modulation/Demodulation Library. 



13 





The Coherent MFSK Demod block demodulates the input, 
which is a complex analog signal (because of the baseband 
implementation of the passband modulation) and is defined by 
five parameters : [M-ary Number, Tone Space, Symbol Interval, 
Initial Phase, Sample Time] . The MFSK Baseband Mod block 
accepts a scalar input in the range [0 M-l] and is defined 
by four parameters: [Tone space, symbol interval, initial 
phase, sample time] . 

From the channel sub-library the AWGN channel block is 
selected and copied into the new window. This block adds 
AWGN to the signal and is defined by three parameters. The 
first two specify the mean and the variance of the noise 
output and the third initializes the seed for the Gaussian 
random number generator. Finally, the error-rate meter block 
is copied from the sink sub-library. This block, shown in 
Figure 8, performs the symbol-to-symbol comparison between 
the sender and the receiver, counts symbol errors, and 
computes the symbol error rate. As shown in Figure 8, 50 
symbols have been transmitted and the receiver has detected 
13 symbols incorrectly. So the symbol error-rate is 
13 

— = 0.270833. All the copied blocks are connected in order to 
50 

complete the MFSK system shown in Figure 9. 
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'i is kj: ag , ;^v^c s 




Figure 9. A Simple Model of a Coherent MFSK System. 
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F. MODEL ANALYSIS 

Comunication system simulation can take place upon the 
completion of the communication system model. There are 
basically three different ways of implementing and 
controlling a simulation. The first is to run the simulation 
interactively from the menu bar and observe the signals at 
various inputs and outputs using oscilloscopes while the 
simulation is running. This method is easy to learn and 
simple to use, but it also incurs the highest overhead and 
is therefore the slowest to execute. The second option is to 
run the simulation by executing a .m file. The main 
advantage of this method is that it provides greater 
flexibility than the first method. For example, parameters 
in the blocks can be changed for repeated simulations and 
initial conditions for various blocks can be overridden. 
Running a simulation from a .m file allows one to change one 
parameter (or parameters) each time the simulation is 
repeated and the results of the simulation can be brought 
into Matlab workspace for further analysis. Finally, the 
third possibility of implementing a simulation (which is the 
most complex and flexibile) is based on the fact that every 
Simulink block has a Matlab equivalent as an S-function. 
Each such S-function has the same name as the corresponding 
model and can be called in a variety of ways to provide 
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information about the model it represents. The S-f unctions 
can provide information about the number of inputs, outputs, 
and states (both continuous and discrete) of the model as 
well as the derivatives of the outputs. All of the analysis 
tools provided with Simulink interact with models through S- 
functions. The linmod and trim .m files are examples, which 
use these functions. These are shown in Table 1. 





Modeling and trimming 


Descreption 


1 


linmod 


Generate linearized model of a 
continuous system 


2 


dlinmod 


Generate linearized model of a 
with discrete elements. 


3 


trim 


Trim systems about equilibrium 
point 



Table 1. Partial List of Analysis Functions (from Ref [9]). 

Generally speaking, the three ways to implement/run a 
simulation are not strictly delineated and different usage 
may be appropriate at different stages of model development 
[Ref 2,6, 9] . 

Finally, Simulink provides a number of methods for the 
numerical solution of ordinary differential equations that 
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can be selected prior to simulation and changed even while 
the simulation is running. These numerical methods for 
solving "initial value problems" are shown in Table 2. 





Integration 


Description 


1 


linsim 


Linear systems method 


2 


rk23 


2 na order Runge-Kutta method 


3 


rk45 


4 tn order Runge-Kutta method 


4 


euler 


Euler's method 


5 


adams 


Adams predictor-corrector method 


6 


gear 


Gear' s predictor-corrector method 
for stiff systems 



Table 2. Functions used for Simulation of Ordinary Dif- 
ferential Equations (from Ref. [9] ) . 
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III. MODEL OF A MFSK PASSBAND DIGITAL COMMUNICATION SYSTEM 

AND ITS BASEBAND EQUIVALENT 

A. BASEBAND EQUIVALENCE FOR PASSBAND SIGNALS AND SYSTEMS 

The most convenient way of representing a narrowband 
system for simulation purposes is to work with equivalent 
baseband quantities. Consider a communication system 
operating with 1 GHz and a bandwidth of 20 MHz. According to 
the Nyquist Sampling Theorem, we need a sampling rate of at 
least 2xf max = 2x1,010 MHz to simulate this system. A 
sampling rate of this magnitude is not practical. Since the 
bandwidth of the baseband equivalent is several orders of 
magnitude smaller than the passband signal's center 
frequency (20 MHz versus 1 GHz), we can simulate the 
passband system by establishing its baseband equivalent 
system and simulating the equivalent system using a 20 
Megasamples/second sampling rate instead of the over 2 
Gigasamples/second sampling rate required for the passband 
signal. Therefore, simulating the baseband equivalent of a 
passband signal substantially reduces the required computer 
storage (memory and disk space) and the calculation time. 

A (real) passband signal can be expressed in polar 
(amplitude and phase) form as: 
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x(t) = a(t) cos[27if c t + cp ( t ) ] => 
x ( t ) = Re{a (t) e j<p(t) e j2,tf c t } (3.1a) 

where a(t) is the envelope, f c is the carrier frequency. Re 
is the real operator and (p(t) is the phase of x(t). 

A (real) passband signal can be also expressed in 
Cartesian form: 

x(t) = Xi(t) cos(27if c t) -x e (t) sin(2;if c t) (3.1b) 

where Xi(t) is the in-phase component and x c (t) is the 
quandrature component of x(t). 

The corresponding complex baseband (or lowpass) signal 
is defined as: 



Xlp (t) = a(t) e j<p(t> (3.2a) 

or: 

x LP (t) = xi(t) + j Xq ( t ) (3.2b) 

In order to transform the real passband signal of 
equation (3.1a) to the complex baseband signal of equation 
(3.2b), we suppress the carrier by dropping e j2,tfct . 

The term a(t)e 3<p(t> is called the complex envelope of 
x(t). The baseband nature of x LP (t) and the low sampling 
rate required to satisfy Nyquist's Theorem for narrowband 
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signals is evident. Complex envelope analysis is very useful 
in analyzing the responses of passband systems to passband 
signals. Finding the response of passband systems to 
passband signals is one of the most essential tasks in 
communications system analysis. 

If we denote the Fourier Transform of the passband 
signal x(t) as X(a>) co=2Ttf and the Fourier Transform of its 
complex envelope x c (t) as X c (cd) then: 

X (to) =0.5 X c (to-© c ) + 0.5 X c *(-o-a> c ) (3.3a) 

where * represents the complex congugate. The first term 
represents complex envelope's baseband spectrum scaled by 
0.5 and shifted to the carrier frequency. The second term 
represents the complex envelope's baseband spectrum scaled 

by 0.5, shifted to carrier frequency, folded about ©=0 axis 
(to negative frequencies) and conjugated. Also, we can 
conclude that the first term represents the spectrum of 

X(to) for positive frequencies while the second term 

represents the spectrum of X(co) for negative frequencies. 

Baseband equivalents of passband systems can be 
determined in a similar manner. For example, consider a 
filter with a passband frequency response H(co) and real 
impulse response h(t). We can now express the frequency 
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response H (co) in terms of the baseband equivalent frequency 
response : 



H (co) =0 . 5 H c (G)-(D C ) + 0.5 H c * (-CO-COe) (3.3b) 

From equations (3.3a) and (3.3b), we can conclude that 
the complex envelope of the passband response to a passband 
input can be obtained as the response of an equivalent 
baseband system to the complex envelope of the passband 
input x(t). In this manner, the baseband equivalent system 
can be used to avoid calculations involving high frequency 
carriers [Ref 5,8]. 

B. MONTE CARLO SIMULATION 

Symbol error probabilities (SEP) and bit error 
probabilities (BER) can be estimated using simulations with 
very large numbers of symbols transmitted by counting the 
symbol errors and taking the ratio of the number of symbol 
errors to the number of symbols as the number of symbols 
tends to infinity, assuming that the symbols are generated 
independently of one another. This is the essence of the 
Monte Carlo Simulation technique. The technique is simple 
but requires significant computing power and/or long 
simulation times when accurate estimates are needed for 
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small symbol/bit error rates. For example, if we count 10 
errors, we can find from Figure 10 the SEP intervals that 
contain the true SEP with 90%, 95%, and 99% confidence [Ref 
5] . 




Figure 10. Confidence bands on P b when observed value is 10 v 
for the Monte Carlo technique (source: Jeruchim etal.,1984). 

From Figure 10, if a symbol (or bit) error probability 
is 10 _1 (v=l), then for 100 experimental errors the actual 
("true") probability of error is between 0.1 to 0.065 and 
0.1 to 0.3 with 95% confidence. However, if errors are not 
independent from each other (errors occur in a "burst") , the 
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first error in a burst is more informative of the error 



statistics than all subsequent errors in the same burst. 
Therefore, in the "dependent error" case we need to count a 
larger number of errors for the same SEP confidence levels, 
which means even longer simulation times [Ref. 5] . If the 

number of experimental errors is denoted r| e and the number 
of the received bits (or symbols) is denoted N b , then the 
confidence with which the BER is upper and lower bounded by 

*i 0 n L 

and , respectively, is given by 

b b 



Prob < BER < — ^-/ti } = © / 3 4) 

N b N b ' 

where @(0<®<1) is the confidence level and r| L , r| D is the 
confidence interval. For a given number of experimental 
errors and if certain reasonable conditions are fulfilled 

(effects of disturbances are independent for each error and 

error-count statistics can be approximated as Gaussian 

distributed) , then the lower and upper number bounds for the 
errors with 95.4% confidence interval are given from the 
following equations: 

r| L =t|e +2-2 yfr\ e +1 (3.5) 
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riu^e +2+2.y/n c +1 



(3.6) 



where r|e is the number of experimental errors and Th,r| 0 is 
the confidence interval. Generally speaking, error bounds 
converge slowly to each other (that is the estimates become 
more accurate) as the number of observed experimental errors 

r| e increases. For ten experimental errors, the true number 
of errors is between 5.4 and 18.6 (between approximately - 
50% and +76% of the observed value), with 95.4% confidence, 
while for 100 experimental errors the true results will be 
between 81.9 and 122 (between approximately -18% and +22% of 
the observed value) also with 95.4% confidence. The values 
of % and r|o from (3.5) and (3.6), respectively, and their 
ratio are tabulated in Table 3 for several values of r| e and 
the confidence level © = 0.954 (95.4% confidence) [Ref. 6]. 
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Experimental 

errors 


Lower limit 
of true er- 
rors (t| L ) 


Upper limit of 
true errors (t| 0 ) 


Ratio-!!b 

% 


1 . 


5 


2.1 


11.9 


5.7 


2. 


10 


5.4 


18.6 


3.5 


3. 


20 


12.8 


31.2 


2.4 


4 . 


50 


37.7 


66.3 


1.8 


5. 


100 


81.9 


122 


1.5 


6 . 


200 


174 


230 


1.3 



Table 3. 95.4% Confidence Intervals for Various Number of 
Experimental Errors (from Ref. [6]). 



C. M-ARY FREQUENCY SHIFT-KEYING (MFSK) 

M-ary frequency-shift keying (MFSK) modulation does not 
require phase coherent detection. This means that MFSK can 
be used in applications where the phase shift incurred by 
the transmission channel varies rapidly, as is the case, for 
example, in communications with high-speed aircraft. With 



MFSK modulation, 


each 


symbol is 


assigned a 


specific 


frequency and an 


MFSK 


signal can 


be described 


by the 



following equation: 



s MFSK (t) =£ A p(t-nT) cos [2Jt(f 0 +a n Af) +cp n ] (3.7) 
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where : 



Af is the spacing between any two adjacent frequencies: 



a n is determined by the information sequence and takes M 

integer values 

f 0 is the lowest frequency 

A is the amplitude of the MFSK signal 

cp n is the phase angle 

p(t-nT) is the unit-amplitude rectangular pulse. 

1 . Bit Error Probability For Coherent MFSK 

In order to obtain the exact expression for the MFSK 
bit error probability, the exact expression for the MFSK 
symbol error probability must be first derived. The coherent 
MFSK symbol error probability is given by the following 
expression [Ref 7]: 



Af=f n -f n _! 



(3.7a) 




(3.8) 



where : 
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£=z/a 



(3.8a) 



a 2 is the variance of wideband Gaussian noise 
M is the number of symbols 
SNR sym is the symbol SNR 



log(M) 

SNR sym = SNRbit7~~ 
log(2) 



( 



and Q(Q is the Q-function given by the following: 



Q (x) = 1-erf {-j= ) ) 



(3 



where the error function (erf) is defined as 



erf(z) = 



1 Z 2 
Je" x dx 



o 



(3 



The following auxiliary functions are introduced 
convenience : 

r(M,SNR Mt ) - Ji^SNR. < 

1 -(C-Y(M,SNR bh )) 2 

v|/(;,M,SNR bit ) = -j= (l-Q(C) ) M_1 e 2 < 3 

V2^ 



.8b) 



. 8c) 



. 8d) 



for 



3.9) 



. 10 ) 
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By using these auxiliary functions we now can express the 
bit error probability [Ref 7] as a function of the number of 
symbols M and the bit signal-to-noise ratio SNR bit : 



Pe (bit) 



1 M r ,SNR blt)+5 

2 M-l Jy(M ; SNR blt )-5 



M/(;,M,SNR blt )dC 



(3.11) 



2 . Bit Error Probability For Noncoherent MFSK 

In order to obtain the exact expression for the MFSK 
bit error probability, the exact expression for the MFSK 
symbol error probability must be first derived. The symbol 
error probability for noncoherent MFSK is given by the 
following expression [Ref 7]: 



• e (sym) 



M-l 

- i-i- 



(— 1) K (M — 1)! 

£o (k + l)k!(M - 1 - k)! ' 



k+12.N„ 



(3.12) 



where : 

M is the number of symbols, M = 2 k , k is an integer, 

N n 2 

—£• is a two-sided power spectral density, N 0 = o 4T, 



and 



A is the amplitude of the MFSK signal. 

We can also express the noncoherent MFSK symbol error 
probability in terms of the symbol energy: 
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(3.13) 



M-l 



■ e (sym) 



= 1-Z 



(-If (M-l)! 

(k + l)k!(M — 1 — k)! 



~k , ^sym 

-e k+1 N ° 



or in terms of the symbol signal-to-noise ratio: 



P 



e (sym) 



! y (-lf(M-l)! ,^snr, 
5J(k + l)k!(M-l-k)! * 



(3.14) 



The symbol error probability for noncoherent MFSK is 
related to the bit-error probability for M = 2 k by [Ref 7] : 



M-l 



1 M n ^ 

Pe(bit) (M,SNR bit ) = -7— -[l-X 

• — 1 k=0 



2 M 



(-lf(M-l)! 



(k + l)k!(M — 1 — k) ! 



] (3.15) 



or 



Pe (bit) ( M, SNRbit ) 



I M ^ (~1) k (M - 1)! L n ( S )SNI " , 

2M-1 fcf(k + l)k!(M-l-k)! 



(3.16) 



D. MFSK MODELS AMD BLOCK ANALYSIS 

1 . Coherent MFSK Model (Noise Only) 

The following figure shows the coherent MFSK (noise 
case only) model, which has been used for the simulation. 
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Figure 11. Coherent MFSK Model for Noise Only Case. 

2 . Block Analysis 

a) Coherent: MFSK Mod Baseband 

Four parameters define the coherent MFSK Mod 
Baseband block: 

[Tone space, symbol interval, initial phase, sample time] . 
The MFSK Baseband Mod block accepts a scalar input in the 
range [0, M-l] . (For example, M in the coherent BFSK model 
is 2) . The output of the MFSK Mod block is a unit-amplitude 
FSK modulated analog signal. The following must be specified 
for the MFSK Mod block: 
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• Tone space (A f ) (also known as the frequency 
separation between two neighboring frequencies and 
defined as A f =l/(2Ts), where T s is the symbol 
interval) . 

• Symbol interval T s (also the data sampling time) 

• Initial phase (the initial phase of the carrier, 
commonly set to 0) 

• Sample frequency (the inverse of the signal sampling 
interval) 

The theoretical BER for coherent BFSK is given as: 



where E B is the average bit energy and N 0 is the single- 
sided noise power spectral density. 

From the Nyquist Theorem, at least two samples per 
shortest period are needed in order to represent a signal. 
For the baseband equivalent model of binary FSK (BFSK) , one 
symbol is represented as DC and the other symbol as a 
carrier with the frequency Af. We have selected the symbol 
interval T s as 1 second. (In the binary case this is also 
the bit interval T b ) . For a good representation of the 
signal we sample at twice the Nyquist rate; that is, four 




(3.17) 
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samples per bit interval, or four samples per second, are 
provided. The number of samples per second N s is given as: 




(3.18) 



where A t is the MFSK signal (not data) sampling time. For 
BFSK and sampling at twice the Nyquist rate, we get: 



For orthogonal MFSK and coherent detection, the frequency 
separation must be an integer multiple of T s /2. The lowest 
separation for orthogonality with coherent detection is: 



This value is defined as the "tone space" (in Hz) in the 
Simulink model. In the case of a coherent 4FSK or 8FSK 
signal the sample time At must be changed because the 
highest frequency in the spectrum is approximately (M-l)Af 
since the maximum frequency for MFSK is defined as f max = (M- 
1) A f . In order to have a good representation of the signal. 




(3.19) 



or 



A t = lsec/4samples per sec A t = 0.25Hz 




=> A f = 0.5Hz 



(3.20) 
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we transmit four samples per symbol interval (or four 
samples per second since we have selected T s =l second) , and 
the sampling frequency for MFSK is: 



f sampling ~ 4 (M 1 ) Af 



(3.21) 



The sampling time At is the reciprocal of the sampling 
frequency: 



At 



sampling 



= l/f 



sampling 



(3.22) 



b) Coherent MFSK Dem Baseband 

This block demodulates the input, which is a 
modulated complex analog signal. Its output is an integer in 
the range [0,M-1]. Five parameters define the Coherent MFSK 
Dem Baseband block: 

• [M-ary number, tone space, symbol interval, initial 
phase, sample time] 

The above parameters should match the ones used in the 
corresponding Coherent MFSK Mod Baseband block. We have 
selected the following values: 



• Tone space (A f ) 



A f = — T Sf where T s is the symbol interval 



Symbol interval (T s ) : 1 sec. 
Initial phase: 0 
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• Sample time : l/(4(M-l)Af) 

• M-ary number: 2, 4, or 8 

c) AWGN Channel 

This block adds AWGN to the signal being 

transmitted through this channel. Three parameters define 
the AWGN channel. The first two specify the mean and the 
variance of the noise output. The third one initializes the 
seed. The seed effectively selects a pseudorandom sequence 
to be generated (for the same value of the seed the same 
sequence is always generated) . 

d) Symbol Counter 

The counter block counts the number of symbols. 

There are three input ports to the counter. The first port 

inputs a constant value of 1. The second port inputs the 

clock signal, which has one pulse per each symbol generated, 
and the third input port is unused ("grounded") . The value 
of 0.5 is set as threshold. Each time the clock pulse 
crosses the threshold going from 0 to 1 (on the rising edge 
of the clock pulse) the counter increments its output by 1, 
thus effectively counting the input symbols. 
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e) Error Counter 



This block compares the delayed "original" signal 
(first input) and the demodulated signal (second input) at 
the end of each symbol interval. The timing is controlled 
by the third input when the rising edge of the square wave 
(from the pulse generator) crosses the threshold. The fourth 
input is not used ("grounded") . The output of this block 
increases by one each time the absolute value of the 
difference between the inputs 1 and 2 exceeds the specified 
tolerance, that is each time the detected symbol differs 
from the original (sent) symbol. Therefore, by setting the 
tolerance of the input difference to 0 (or a very small non- 
zero value) , the symbol errors are counted. 

f ) Triggered Write To Pile 

Four parameters define the Triggered Write to File 

block: 

[File name, data type, number of trigger pulses between 
saved data, threshold in detecting trigger signal] 

This block writes a record to a file only at the rising edge 
of the signal coming from the pulse generator (the block can 
also isolate and write all the errors to a file in case of 
burst errors) . The file can be an ASCII text file or a data 
type ("integer," "float"). The third parameter (number of 
trigger pulses between saved data) determines the number of 
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rising edges that should be received between file writes. 
Setting this field to zero forces the block to write a 
record at every rising edge of the trigger. This block has 
two input ports. The first input is the error counter value 
that is recorded into a file and the second input input is 
the "write trigger" signal from the pulse generator. 

g) Triggered Write To Workspace: 

Six parameters define the Triggered Write to 
Workspace block: 

[Workspace variable, data type, number of trigger pulses 
between saved data, maximum row number, section to keep if 
overflow occurs, threshold in detecting trigger signal] 

The saved workspace variable can be either a column vector 
when the first input port (message signal) has a scalar 
input or a matrix when the first input has a vector input. 
The first element of the input signal vector is saved in the 
first column; the second element of the input signal vector 
is saved in the second column, and so on. The block can save 
the input as string variables or data. The maximum row of 
the output variable is limited by a pre-defined number, 
which can be changed during simulation. After the limit is 
reached, the block keeps the portion of the data defined by 
the "overflow" entry. Data may be written into the workspace 
at every trigger or a number of triggers may be skipped 
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between sucessive writes, depending on the value of the 
corresponding entry (zero forces a write at every trigger) . 

h) K-Step Delay 

This block implements a delay that is equal to an 
integer multiple of the sampling interval that is in effect 
for the particular simulation. 

i) Sum 

This block adds or subtracts the two inputs 
depending on the sign selection for the inputs: ++, +-, -+, 
Since the objective is to compare the original and the 
detected symbols, either +- or -+ can be used. 

j) Abs 

This block's output is the absolute value of the 
input. This block creates a positive value of 1 each time 
the original and the detected symbols differ, regardless of 
the sign (polarity) of the difference. 

k) Pulse Generator 

Four parameters define the pulse generator block: 
[Pulse period (symbol period), duty cycle (% of period), 
amplitude, start time] 

This block generates a unit-amplitude square wave with the 
period equal to the symbol interval. The pulse generator 
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starts at the end of the first symbol interval such that the 
rising edges of the square wave coincide with the end of 
symbol times. Therefore, since there is a rising edge of the 
square wave at the end of each symbol interval, there is 
always a rising edge coincident with a symbol decision time. 

l) Switch 

The switch block has three input ports. The first 
port inputs the clock signal and the second input port 
inputs the absolute value of the sum of the delayed and 
demodulated signal. The third input is unused ("grounded"). 

m) Sampled Read From Workspace 

The Sampled ' Read From Workspace block reads a 
matrix from MATLAB workspace. In our case the matrix is a 
matrix of pseudo-random integers generated using MATLAB' s 
random integer function randint. There are four argument 
fields, delineated by commas, in the specification of the 
randint function. The first two define the size (rows and 
columns) of the matrix to be generated. The third field 
accepts either a single integer (the number of symbols M) or 
two-integer elements array [0, M-l] . In either case the 
range of random integers (the elements of the random matrix) 
is from 0 to M-l and all M values are equally probable. The 



39 



last field is the seed for the pseudorandom matrix (for the 

same value of the seed the same matrix is always repeated) . 

Randint [number of rows, number of columns, [0 number of 
symbols -1], seed] 

or 

Randint [number of rows, number of columns, number of 

symbols, seed] 

In the case of noncoherent MFSK (with noise only) the 
same model as for the coherent MFSK (noise case only) is 
used with only one difference. The coherent MFSK Mod 

Baseband and coherent MFSK Dem Baseband blocks are replaced 
with the noncoherent MFSK Mod Baseband and noncoherent MFSK 
Dem Baseband blocks. These blocks are defined by the same 
parameters as in the coherent MFSK Mod Baseband and coherent 
MFSK Dem Baseband blocks. Also, the frequency separation is 

Af = 1/T S instead of 1/ (2T S ) . 

3 . Coherent MFSK Model For MFSK With Interference And 
Additive Noise 

The coherent MFSK model with interference AWGN is shown 
in Figure 12. 
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Figure 12. Model for Coherent MFSK with Interference and 
Additive Gaussian Noise. 

The difference between this model and the model for 
coherent MFSK with noise only is that the white blocks 
representing the interference signal have been added. The 
interference signal is assumed synchronized with the desired 
signal. That is, symbol transitions of the interference 
occur at the same instants in time as the symbol transitions 
of the desired signal. For noncoherent MFSK with 
interference and AWGN, the same model is used (Figure 13) 
with the coherent MFSK Mod Baseband and coherent MFSK Dem 
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Figure 13. Model for Noncoherent MFSK with Interference 
and Additive Gaussian Noise. 

Baseband blocks replaced with the noncoherent MFSK Mod 
Baseband and noncoherent MFSK Dem Baseband blocks and the 
frequency separation increased to 1/T S . In order to run 
multiple simulations for a number of different values of 
parameters such as the signal-to-noise ratio (SNR) , Matlab 
.m files were used. For coherent MFSK, the Matlab programs 
sun_prepl.m and sun_fskl.m (listed in the Appendix) were 
used. These programs can also be combined into a single 
program. For noncoherent MFSK, essentially the same Matlab 
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.m files were used, with the appropriate variables changed. 
The .m files are used to complete the three basic tasks: 
data generation and input, simulation, and recording the 
output into files, as shown in Figure 14. 




Figure 14. Flow Chart for Matlab .m File Simulation 
Implementation . 

The first task involves the choices concerning the case 
we want to test (noise or noise and interference), the type 
of MFSK (2FSK, 4FSK, 8FSK...) , and the simulation parameteres 
such as symbol duration T s , oversampling factor, minimum and 
maximum signal-to-noise and signal-to-interf erence ratios, 
mimimum number of errors acceptable, the factor multiplying 
the number of errors (to select the number of symbols to 
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try) , the maximum size of the random integer array, and the 
name of the file for saved data. The second task is the 
actual simulation started by the appropriate Matlab command 
and using the selected Simulink model (coherent MFSK or 
noncoherent MFSK) . The third task is recording of the 
results: the signal-to-noise and signal-to-interference 

ratios, the number of errors, the number of symbols, and 
symbol error probability, all as ASCII data files. 
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IV. SIMULATION, ANALYSIS, AND PERFORMANCE VERIFICATION FOR 

AWGN 

A. PERFORMANCE VERIFICATION FOR AWGN 

In this chapter computer simulation results are 
presented in order to verify the performance of MFSK 
communication systems in the presence of AWGN. These results 
are the average of three computer simulation runs for each 



signal- 


to-noise (SNR) 


ratio . 


The value 


of two 


has 


been 


used 


as an 


oversampling 


factor 


(sampling 


rate 


of 


twice 


the 


minimum 


value) . The 


simulations ran 


until 


at 


least 


100 



errors were observed. The bit error probability intervals 
for the 90%, 95%, and 99% confidence levels may be found 

from Figure 10 (Chapter III) . For a given number of 

experimental errors and for the confidence level of 95.4%, 
the lower and upper bounds on the number of errors can be 
found from Table 3-1. The observed error number is a random 
variable, subjected to statistical scatter arising from 
different realizations (in general, a different number of 
errors is observed each time a simulation is run with a 
different set of seeds for data and noise) . In order to 

prevent 'out of memory' errors, the data sequences were 

limited to 10 6 symbols for each simulation. If less than 100 
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errors were observed in a simulation, then the sequence was 
repeated until a sufficient number of errors occurred. 



1 . Results For Coherent Binary FSK (BFSK) 

The theoretical and the experimental bit error 
probability for coherently detected binary FSK (BFSK) as a 
function of the SNR in dB is shown in Figure 15. 
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Figure 15. Probability of Bit Error (Theory and Simulation) 
for Coherent BFSK Detection. 
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The relative difference between the theoretical 



probability of error for coherent BFSK and the probability 
of error estimates obtained by the simulation are shown in 
Figure 16. 
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Figure 16. Percent difference between Theoretical and 
Simulation Results for Coherent Detection of BFSK. 

The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 



Percent Error 



! ; ?! 

t ' 

» * i ! 

• \ 

• • 1 

Q. 


'• I 

*. ! 

0 — o 

# • » 


1 _ 

7 

| 




/ • ! ; 
' 1 * / 

* • t 

O i 

/ * • 

• • , * 




# ; •' • 
i | i « i 

# 1 i * * 

gl. / ■ r 

$ * M 

• -<3 © 

. __ # ..... t ...... .1 . . . ........ ... 


! 


1 

» 

* 


1 « 
* 


i r 1 ; 

; i 1 

! | ! 

i ; | i 

! i j : 

! ! j 


I 


• 

1 

• 

» 

© 


* : : i ; j 

i 1 : ; 

• 

j • j i • i i 

1 i | ; 1 





47 



Figure 17. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that minimum has been reached. For this particular 
simulation, convergence is achieved at approximately 1.2'10 5 
symbols . 
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Figure 17. Convergence Plot for Coherent Detection of BFSK. 
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The accuracy of the simulation estimate for the bit error 
probability for coherent BFSK is shown in Figure 18. The 
convergence to zero as the number of symbols increases is 
evident . 




P,„ 2»10 4 4»10 4 6-1 0 4 8-10 4 1-10 5 1 .2-1 0 5 1.4-10 5 , 

Sym 2 n J.25865-10 5 , 

Figure 18. Percent Difference Between Theory and Simulation 
for Coherent BFSK as a function of the Number of Symbols. 

The differences between the theoretical bit error 
probabilities for coherent 2FSK and the simulation estimates 
for various SNR's are plotted in Figure 19 as a histogram. 
All estimates are within the -20% to +10% interval about the 



49 



theoretical values, with most differences very close to 0 
(approximately +2% off) . 
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Figure 19. Histogram of Simulation Errors for Coherent BFSK. 

2 . Results For Coherent 4FSK 

The theoretical and the experimental bit error 
probability for coherently detected 4FSK as a function of 
the SNR in dB is shown in Figure 20. 
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Figure 20. Probability of Bit Error (Theory and Simulation) 
for Coherent 4FSK Detection. 

The relative difference between the theoretical 
probability of error for coherent 4FSK and the probability 
of error estimates obtained by the simulation is shown in 
Figure 21. 
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Figure 21. Percent difference between Theoretical and 

Simulation Results for Coherent Detection of 4FSK. 

The convergence of the estimate of the bit error 
probability as a function of the number of symbols 

transmitted in the course of the simulation is shown in 
Figure 22. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
simulation, the convergence is evident at approximately 

4.5‘10 4 symbols. 
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Convergence 




Figure 22. Convergence Plot for Coherent 4FSK. 

The accuracy of the simulation estimate for the bit error 
probability for coherent 4FSK is shown in Figure 23. The 
'oscillatory' convergence to zero as the number of symbols 
increases is evident. 
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Figure 23. Percent Difference Between Theory and 
Simulation for Coherent 4FSK as a function of the 
Number of Symbols. 

The differences between the theoretical bit error 
probabilities for coherent 4FSK and the simulation estimates 
for various SNRs are plotted in Figure 24 as a histogram. 
All estimates are within the -12% to +40% interval about the 
theoretical values, with most differences between +10% and 
+30%. 
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Figure 24. Histogram of Simulation Errors for Coherent 
4FSK. 

3 . Results For Coherent 8FSK 

The theoretical and the experimental bit error 
probability for coherently detected 8FSK (8FSK) as a 
function of the SNR in dB is shown in Figure 25. 
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Figure 25. Probability of Bit Error (Theory and 
Simulation) for Coherent 8FSK Detection. 

The relative difference between the theoretical 
probability of error for coherent 8FSK and the probability 
of error estimates obtained by the simulation is shown in 
Figure 26. 
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Figure 26. Percent difference between Theoretical and 
Simulation Results for Coherent Detection of 8FSK. 

The convergence of the estimate the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 27. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
simulation, the convergence is evident at approximately 
2.5'10 4 symbols. 
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Convergence 




Figure 27. Convergence Plot for Coherent Detection of 8FSK. 

The accuracy of the simulation estimate for the bit 
error probability for coherent 8FSK is shown in Figure 28. 
The 'oscillatory' convergence to approximately zero as the 
number of symbols increases is evident . 
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Figure 28. Percent Difference Between Theory and 
Simulation for Coherent 8FSK as a function of the Number 
of Symbols. 

The differences between the theoretical bit error 
probabilities for coherent 8FSK and the simulation estimates 
for various SNRs are plotted in Figure 29 as a histogram. 
All estimates are within the -22% to +0% interval about the 
theoretical values, with most differences between -7% and 

_ o_ 



59 



0.25 



Histogram 



0.2 n ; irn t 

0 .15 : ; : ; ! 

hist error 
4- " q 

0.1 ; ; 1 n fl 

0.05 1 : j ? 

A 

0 1 11 L “ 1 1 11 

-25 -20 -15 "10 "5 0 5 

2 1.892 182 centers 5 

Figure 29. Histogram of Simulation Errors for Coherent 
8 FSK . 

4 . Results For Noncoherent BFSK 

The theoretical and the experimental bit error 
probability for noncoherently detected binary FSK (BFSK or 
2FSK) as a function of the SNR in dB is shown in Figure 30. 
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Figure 30. Probability of Bit Error (Theory and 
Simulation) for Noncoherent BFSK Detection. 



The relative difference between the theoretical 
probability of error for noncoherent BFSK and the 
probability of error estimates obtained by the simulation is 
shown in Figure 31. 
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Figure 31. Percent difference between Theoretical and 
Simulation Results for Noncoherent Detection of BFSK. 



The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 32. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
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simulation, the convergence is evident at approximately 
2.3'10 4 symbols. 




Figure 32. Convergence Plot (100 errors) for Noncoherent 
Detection of BFSK. 

The accuracy of the simulation estimate for the bit 
error probability for noncoherent BFSK is shown in Figure 
33. In this case, the convergence is to a non-zero value 
(approximately -30% off) as the number of symbols increases. 
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Figure 33. Percent Difference Between Theory and Sim- 
ulation for Noncoherent BFSK as a function of the Num- 
ber of Symbols. 



The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation using 200 
minimum acceptable errors is shown in Figure 34. This allows 
us to estimate the minimum number of symbols required for a 
particular simulation and to verify that the minimum has 
been reached. For this particular simulation, the 
convergence is evident at approximately 8'10 4 symbols. 
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Figure 34. Convergence Plot (200 errors) for Noncoherent 
Detection of BFSK. 

The accuracy of the simulation estimate for the bit 
error probability for noncoherent BFSK is shown in Figure 
35. In this case, the convergence is again to a non-zero 
value (approximately -30% off) as the number of symbols 
increases . 
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Figure 35. Percent Difference Between Theory and Simulation 
for Noncoherent BFSK (200 Errors) as a function of the 
Number of Symbols. 

The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation using 300 
minimum acceptable errors is shown in Figure 36. This allows 
us to estimate the minimum number of symbols required for a 
particular simulation, and to verify that the minimum has 
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been reached. For this particular simulation 
convergence is evident at approximately 9‘10 4 symbols. 



the 




Figure 36. Convergence Plot (300 errors) for Noncoherent 
Detection of BFSK. 



The accuracy of the simulation estimate for the bit 
error probability for noncoherent BFSK is shown in Figure 
37. In this case, the convergence is again to a non-zero 
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value (approximately -30% off) as the number of symbols 



increases . 




Figure 37. Percent Difference Between Theory and 
Simulation for Noncoherent BFSK (300 Errors) as a function 
of the Number of Symbols. 

The differences between the theoretical bit error 
probabilities for noncoherent BFSK (100 errors case) and the 
simulation estimates for various SNRs are plotted in Figure 
38 as a histogram. All estimates are within the -22% to +87% 
interval about the theoretical values, with most differences 
between -22% and 0%. 
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Figure 38. Histogram of Simulation Errors (100) for 
Noncoherent BFSK Detection. 



5. Results For Noncoherent 4FSK 

The theoretical and the experimental bit error 
probability for noncoherently detected 4FSK as a function of 
the SNR in dB is shown in Figure 39. 
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Figure 39. Probability of Bit Error (Theory and 
Simulation) for Noncoherent 4FSK Detection. 

The relative difference between the theoretical 
probability of error for noncoherent 4FSK and the 
probability of error estimates obtained by the simulation is 
shown in Figure 40. 
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Figure 40. Percent difference between Theoretical and 
Simulation Results for Noncoherent Detection of 4FSK. 

The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 41. This allows us to estimate the minimum number of 
symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
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simulation, the convergence is evident at approximately 
9'10 3 symbols. 




Figure 41. Convergence Plot for Noncoherent Detection of 
4 FSK. 

The accuracy of the simulation estimate for the bit 
error probability for noncoherent 4 FSK is shown in Figure 
42. In this case the convergence is again to a non-zero 
value (approximately -20% off) as the number of symbols 
increases . 
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Figure 42. Percent Difference Between Theory and 
Simulation for Noncoherent 4FSK as a function of the 
Number of Symbols. 

The differences between the theoretical bit error 
probabilities for coherent 4FSK and the simulation estimates 
for various SNRs are plotted in Figure 43 as a histogram. 
All estimates are within the 5% to 70% interval about the 
theoretical values, with most differences between 5% and 
20 %. 
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Figure 43. Histogram of Simulation Errors for Noncoherent 
4FSK. 



6 . Results For Noncoherent 8FSK 

The theoretical and the experimental bit error 
probability for noncoherently detected 8FSK as a function of 
the SNR in dB is shown in Figure 44. 
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Figure 44. Probability of Bit Error (Theory and Simula- 
tion) for Noncoherent 8FSK Detection. 

The relative difference between the theoretical 
probability of error for noncoherent 8FSK and the 
probability of error estimates obtained by the simulation is 
shown in Figure 45. 
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Figure 45. Percent difference between Theoretical and 
Simulation Results for Noncoherent Detection of 8FSK. 

The convergence of the estimate of the bit error 
probability as a function of the number of symbols 
transmitted in the course of the simulation is shown in 
Figure 46. This allows us to estimate the minimum number of 
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symbols required for a particular simulation and to verify 
that the minimum has been reached. For this particular 
simulation, the convergence is evident for approximately 
5'10 3 symbols. 




Figure 46. Convergence Plot for Noncoherent Detection of 
8FSK. 

The accuracy of the simulation estimate for the bit error 
probability for noncoherent 8FSK is shown in Figure 47. In 
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this case, the convergence is again to a non-zero value 
(approximately -13% off) as the number of symbols increases. 
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Figure 47. Percent Difference Between Theory and 
Simulation for Noncoherent 8FSK as a function of the 
Number of Symbols. 



The differences between the theoretical bit error 
probabilities for coherent 4FSK and the simulation estimates 
for various SNR's are plotted in Figure 48 as a histogram. 
All estimates are within the -12% to +5% interval about the 
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theoretical values, with most differences between -7% and 
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Figure 48. Histogram of Simulation Errors for 
Noncoherent 8FSK. 



B . OBSERVATIONS 

Simulation results for coherent MFSK in the presence of 
AWGN agree very well with theory. As can be seen from the 
convergence plots, in all cases of coherent MFSK detection 
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the differences between the bit error probability estimates 
and the theoretical values converge to zero as the number of 
symbols increases, which means that the simulation estimates 
are correct. However, this is not the case for noncoherent 
MFSK . The simulation results are too small for BFSK and 8FSK 
(bit error probability is underestimated) compared to the 
theoretical ones. Only in case of 4 FSK are the simulation 
results higher (bit error probability is overestimated) 
than the theoretical ones. The simulation results converge 
to values that are approximately: 

• 30% lower for noncoherent BFSK 

• 20% higher for noncoherent 4FSK and 

• 13% lower for noncoherent 8FSK. 

The simulation results for both noncoherent 4FSK and 
8FSK are within expected bounds. Only noncoherent BFSK is 
not within expected bounds for the 95.4% confidence 

intervals for 100 errors of Monte Carlo Simulations. An 
increase in the observed number of errors to 200 and 300 
does not alter the simulation results, which still converge 
to the same somewhat incorrect values. This can be seen in 
Figure 33, 34, 35, and 36 for noncoherent BFSK where, 

despite the fact that 100, 200, and 300 errors were 
observed, the bit error probability still converged to about 
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the same (incorrect) value. This suggests a systematic error 
in the Communications Toolbox for noncoherent BFSK. 



The summary of the accuracy of the simulation for both 
coherent and noncoherent MFSK detection is presented in 
Tables 4,5, and 6. 





Bit Error Probability 


Coherent BFSK 


Noncoherent BFSK 


1 . 


Mean Percent Error 


-2.532 


0.538 


2. 


Max. Percent Error 


9.065 


87.28 


3. 


Min. Percent Error 


-18.209 


-21.12 


4 . 


Standard Deviation 


7.777 


31.571 



Table 4. Comparison of the Simulation and the Theoretical 
Results for Coherent and Noncoherent BFSK Detection. 



From Table 4, we note that the standard deviation is 
much larger for the noncoherent MFSK. This means that the 
simulation results are less accurate for noncoherent BFSK 
detection . 





Bit Error Probability 


Coherent 4FSK 


Noncoherent 4FSK 


1 . 


Mean Percent Error 


13.628 


18.946 


2. 


Max. Percent Error 


38.725 


4.101 


3. 


Min. Percent Error 


-10.883 


69.43 


4. 


Standard Deviation 


14.109 


17.931 



Table 5. Comparison of the Simulation and the Theoretical 
Results for Coherent and Noncoherent 4FSK Detection. 
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From Table 5, we note that the mean values and the 
standard deviations for coherent and noncoherent 4FSK are 
closer to each other (the mean values in the absolute 
sense) , so the accuracy of the simulation is comparable for 
both coherent and noncoherent 4FSK. 





Bit Error Probability 


Coherent 8FSK 


Noncoherent 8FSK 


1 . 


Mean Percent Error 


-8.207 


-1.666 


2. 


Max. Percent Error 


0.191 


21.585 


3. 


Min. Percent Error 


-20.105 


-10.704 


4. 


Standard Deviation 


6.281 


8.873 



Table 6. Comparison of the Simulation and the Theoretical 
Results for Coherent and Noncoherent 8FSK Detection. 



From Table 6, we note that the mean values and the 
standard deviations for coherent and noncoherent 8FSK are 
close to each other (the mean values in the absolute sense) , 
so the accuracy of the simulation is in general comparable 
for both coherent and noncoherent 8FSK. It is interesting to 
note that as the number of symbols M increases from 2 to 8 
the simulation results for noncoherent MFSK become more 
accurate while the simulation results for coherent MFSK 
become less accurate. 
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V. SIMULATION AND PERFORMANCE ANALYSIS FOR CO-CHANNEL 

INTERFERENCE 

A. INTERFERENCE IN DIGITAL COMMUNICATION SYSTEMS 

Interference is a factor limiting the performance of 
digital communication systems. Co-channel interference is 
one of major types of system-generated interference and 
refers to the degradation caused by an interfering waveform 
appearing within the signal bandwidth. It can be introduced 
in a variety of ways, but most commonly by other users of 
the same portion of RF spectrum operating similar types of 
equipment [Ref 3] . In Chapter 3 MFSK digital communication 
systems were modeled, and in Chapter 4 the performance of 
these MFSK systems operating in the presence AWGN was 
verified. Since theoretical results for bit error 
probabilities are available only for the cases of AWGN, one 
must resort to simulations to estimate the bit error 
probabilities for MFSK systems with both AWGN and co-channel 
interference . 

In this chapter computer simulation results are 
provided for BFSK, 4FSK, and 8FSK communication systems 
operating in the presence of wideband, AWGN and co-channel 
interference. The results are the averages of two computer 
simulations for each case. The sampling rate was selected as 
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twice the minimum sampling rate (an oversampling factor of 
two) , and simulations were run until at least 100 errors 
were observed. In order to prevent 'out of memory' errors, 
the data sequences were limited to 10 6 symbols for each 
simulation. If less than 100 errors were observed in a 
simulation, the sequence was repeated until a sufficient 
number of errors occurred. Finally, the computer simulation 
results were obtained for a range of -5dB to +15 dB symbol 
signal-to-noise and symbol signal-to-interf erence (i.e., 
jamming) ratios. 

B. PERFORMANCE OF MFSK DIGITAL COMMUNICATION SYSTEMS IN 
THE PRESENCE OF AWGN AND CO-CHANNEL INTERFERENCE 

In order to specify the performance of MFSK 
communication systems operating in the presence AWGN and co- 
channel interference, two sets of plots are provided for M = 
2, 4, and 8, for both coherent and noncoherent detection. 
The bit error probability as a function of the signal-to- 
noise ratio (SNR) with the signal-to-interf erence ratio 
(SJR) as a parameter is shown in the first set, while the 
second set shows the bit error probability as a function of 
the signal-to-interference ratio (SJR) with the signal-to- 
noise ratio (SNR) as a parameter is shown in the second set. 
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1. 



Results For Coherent Binary FSK (BFSK) 



The probability of bit error for coherent BFSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
49. The BER increase as SNR decreases, which only happens 
for negative SJR, is discussed in Appendix E. 



_0 514167 _ 



e_bit (M 



. SNR vi, 



bit ) 



BER avg m ,„ 

BER avg m , 1 

B-B-B 

BER -g m , 2 

BER 'avg m , 3 

v/v/v/ 

AA A 

BER av* m , 4 

BER -g m>5 

□ D □ 

BER -g m , 6 



avg 



BER 



*■ Iv* 

BER a vg m ,g 

BER avg m>9 

o o o 



,2.660029 10 



-4 



1-10 



1-10 




-5 



0 5 10 

10 -log (SNR bi , j.SNR m 



Figure 49. Probability of Bit Error versus SNR for 
Coherent BFSK with AWGN and Interference. 
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Ten curves for signal-to-interference (SJR) ratios 
from -5 dB to +12 dB in 1.89 dB increments are shown in 
Figure 49. The solid line is the probability of bit error 
for the BFSK without co-channel interference. We note the 
dramatic increase in the bit error probability due to co- 
channel interference relative to the AWGN only. 
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Figure 50. Probability of Bit Error versus SJR for 
Coherent BFSK with AWGN. 
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The probability of bit error for coherent BFSK as a 
function of the signal-to-interference ratio SJR is shown in 
Figure 50. Ten curves for signal-to-noise (SNR) ratios from 
-5 dB to +12 dB in 1.89 dB increments are shown. The solid 
line (unit step function) is the probability of bit error 
for the case when only signal and interference are present 
(AWGN) . Note that although the unit step has the value of 
zero for positive SJR, the value of 10' 4 has been selected 
for SJR > 0 such that a logarithmic scale that shows the 
simulation results well (all > 10" 4 ) can be used. When the 
interference power is less than the signal power, (positive 
SJR) the interference has no effect, and the probability of 
error is zero because the receiver always selects the larger 
of the detected symbols. On the other hand, if the 
interference power is larger than the signal power, the 
probability of error is 0.5 since (statistically) half of 
the interference symbols will be opposite to the signal 
symbols and will be selected by the receiver because of 
their larger power. Consequently, the receiver will be 
making incorrect decisions 50% of the time. 

In the noise-free case, the probability of bit errors 
caused by co-channel interference exhibits step-like 
behavior (from 0 to 0.5) with the threshold at SJR = 0 dB. 
As the signal-to-interference ratio increases (that is, as 
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the interference power decreases relative to the signal 
power) , the probability of bit error decreases at first, but 
then becomes essentially constant (the curves become nearly 
horizontal lines) . This indicates that there are roughly two 
regimes of operation: interference-dominated and noise- 
dominated. Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be defined by 
inspection. For example, the 7-th curve (for SNR=6.333 dB) 
becomes flat for SJR > 12 dB meaning that the performance of 
the system is interference dominated for SJR < 12 dB. 

In general, the condition SNR = SJR seems appropriate 
as the border between the interference-dominated (SNR>SJR) 
and noise-dominated (SJR>SNR) , regimes of operation for 
lower values of the error probability. This is because the 
curves for the bit error probability are steep for higher 
values of SNR and SJR, meaning that a small change in SNR or 
SJR causes a large change in the bit error probability. 
Therefore, for small values of the bit error probability 
(larger values of SNR and SJR) the smaller of the two ratios 
(SNR or SJR) controls the bit error probability. 
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2 . 



Results For Coherent 4FSK 



The probability of bit error for coherent 4FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
51. 
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Figure 51. Probability of Bit Error versus SNR for 
Coherent 4FSK with AWGN and Interference. 



Ten curves for symbol signal-to-interference (SJR) 
ratios from -5 dB to +15 dB in 2.223 dB increments are shown 
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in Figure 51. The solid line is the probability of error for 
4FSK without co-channel interference (AWGN) . We note the 
dramatic increase in the bit error probability due to co- 
channel interference relative to AWGN only. The probability 
of bit error for coherent 4FSK as a function of the signal- 
to-interference ratio SJR is shown in Figure 52. 
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Figure 52. Probability of Bit Error versus SJR for 
Coherent 4FSK with AWGN and Interference. 
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Ten curves for symbol signal-to-noise (SNR) ratios from 
-5 dB to +15 dB in 2.223 dB increments are shown in Figure 
52. The solid line (unit step function) shows the 
probability of bit error for the case when only interference 
is present (no noise) . Since the plots in Figures 50 and 51 
represent the probability of bit error for coherent 4FSK, 
the 'unit step' may also be defined on per-bit basis. For M 
= 4 the unit step is shifted to the left by a factor of 

1 0 • log- ( 1°^) ) w hich is approximately 3 dB. Also, although the 
log(M) 

unit step has the value of zero for positive SJR, the value 
of 1CT 7 is selected such that a logarithmic scale that shows 
the simulation results (which are all > 1CT 7 ) can be used. 

As was found in Section 1 (Results for Coherent BFSK) , 
there are roughly two regimes of operation: interference- 
dominated and noise-dominated. Although the curves do not 
change abruptly from decreasing monotonically to being 
horizontal, approximate values that separate the two regions 
may be again identified by inspection. For example, the 8- 
th curve (for the symbol bit SNR=7.545 dB) becomes flat for 
SJR > 12 dB, meaning that the performance of the system is 
noise dominated for SJR > 12 dB. In general, as found for 
coherent BFSK, SNR = SJR is appropriate as the border 
between the interference-dominated (SNR>SJR) and noise- 
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dominated (SJR>SNR) regimes of operation for low values of 
the bit error probability. Therefore, for low values of bit 
error probability (larger values of SNR and SJR) the smaller 
of the two (SNR or SJR) by and large controls the bit error 
probability. 

3 . Results For Coherent 8FSK 

The probability of bit error for coherent 8FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
53. Ten curves for symbol signal-to-interf erence (SJR) 
ratios from -5 dB to +15 dB are shown. The solid line is the 
probability of error for the coherent 8FSK without co- 
channel interference (AWGN) . We note the dramatic increase 
in bit error probability due to co-channel interference 
relative to AWGN only. 
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Figure 53. Probability of Bit Error versus SNR for 
Coherent 8FSK with AWGN and Interference. 



The bit error for coherent 8FSK as a function of the 
signal-to-interference ratio SJR is shown in Figure 54. 



93 



BER int only \ 1 ® JR 

BERav S0,n 

XXX 

BER avg ! „ 

BBG 

BER .vg 2 ,„ 

ber » 83 ,„ 

vLyLV 
•\ A A 

BER ” S4 , n 



ber ^ 5 ,„ 

BER avg 6>n 

0 - 0 - 0 

BERav S7,n 

v/v/v/ 

A A A 

BERav g g n 

BERav §9,n 

o - e -e 




-5 0 5 10 



-8 



SJR bitdB’ SJR bit_dB n 



13 



Figure 54 . Probability of Bit Error versus SJR for 
Coherent 8FSK with AWGN and Interference. 



Ten curves for symbol signal-to-noise (SNR) ratios from 
-5 dB to +15 dB in 2.223 dB increments are shown in Figure 
54. The solid line (unit step function) is the probability 
of bit error for the case when only signal and interference 
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are present (no noise) . Since the plots in Figures 53 and 54 
represent the probability of bit error for coherent 8FSK, 
the 'unit step' may also be defined on per-bit basis. For M 
= 8 the unit step is shifted to the left by a factor of 

10-loe-(— , which is approximately 4.71 dB. Also, although 

log(M)' 

the unit step has the value of zero for positive SJR, the 
value of 10" 6 has been selected such that a logarithmic 
scale that best shows the simulation results (which are all 
> 10” 6 ) could be used. 

As with the two previous cases, there are roughly two 
regimes of operation: interference-dominated and noise- 

dominated. Although the curves do not change abruptly from 
decreasing monotonically to being horizontal, approximate 
values that separate the two regions may be identified by 
inspection. For example, the 8-th curve (for SNR=5.784 dB) 
becomes flat for bit SJR > 10 dB meaning that the 
performance of the system is noise dominated for bit SJR > 
10 dB. In general, as with the two previous cases, the 
condition SNR = SJR seems to be the border between the 
interference-dominated (SNR>SJR) and noise-dominated 
(SJR>SNR) regimes of operation for lower values of bit error 
probability. Therefore, for small values of the bit error 
probability (larger values of SNR and SJR) the smaller of 
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the two (SNR or SJR) by and large controls the error 
probability. 

4 . Results For Noncoherent Binary FSK (BFSK) 

The probability of bit error for noncoherent BFSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
55. 
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Figure 55. Probability of Bit Error versus SNR for 
Noncoherent BFSK with AWGN and Interference. 
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Ten curves for signal-to-interference (SJR) ratios from 
-5 dB to +12 dB in 1.8 9 dB increments are shown in Figure 
55. The solid line is the probability of bit error for 
noncoherent BFSK without co-channel interference (AWGN 
only) . We note the dramatic increase in the bit error 
probability due to co-channel interference relative to AWGN 
only. 

The probability of bit error for noncoherent BFSK as a 
function of the signal-to-interference ratio SJR is shown in 
Figure 56. Ten curves for signal-to-noise (SNR) ratios from 
-5 dB to +12 dB are shown. The solid line (unit step 
function) is the probability of bit error for the case when 
only interference is present (no noise) . Note that the unit 
step has the value of zero for positive SJR but the value of 
10" 3 has been selected so that a logarithmic scale that 
shows the simulation results (which are all > 10' 3 can be 
used. As with the previous cases examined, there are roughly 
two regimes of operation: interference-dominated and noise- 
dominated. Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be identified by 
inspection. For example, the 7-th curve (for SNR=6.333 dB) 
becomes flat for SJR > 12 dB meaning that the performance of 
the system is noise dominated for SJR > 12 dB . In general. 
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as before, the condition SNR = SJR is the border between the 
interference-dominated (SNR>SJR) and noise-dominated 
(SJR>SNR) regimes of operation for lower values of the error 
probability. 
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Figure 56. Probability of Bit Error versus SJR for 
Noncoherent BFSK with AWGN and Interference. 
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Therefore, for small values of the error probability 
(larger values of SNR and SJR) , the smaller of the two (SNR 
or SJR) by and large controls the error probability. 

5 . Results For Noncherent 4FSK 

The probability of bit error for noncoherent 4FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
57. 
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Figure 57. Probability of Bit Error versus SNR for 
Noncoherent 4FSK with AWGN and Interference. 
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Eleven curves for symbol signal-to-interf erence (SJR) 
ratios from -5 dB to +15 dB in 1.89 dB increments are shown 
in Figure 57. The solid line is the probability of bit error 
for the noncoherent 4FSK without co-channel interference 
(AWGN only) . We note the dramatic increase in the bit error 
probability due to co-channel interference, relative to AWGN 
only. 

The probability of bit error for noncoherent 4FSK as a 
function of the signal-to-interf erence ratio SJR is shown in 
Figure 58. Eleven curves for symbol signal-to-noise (SNR) 
ratios from -5 dB to +15 dB in 1.89 dB are shown. The solid 
line (unit step function) is the probability of bit error 
for the case when only interference is present (no noise) . 
Since the plots in Figures 57 and 58 represent the 
probability of bit error for noncoherent 4FSK, the 'unit 
step' may also be defined on per-bit basis. For M = 4 the 

unit step is shifted to the left by a factor of 1 0 • log- ( \ 

log(M) 

which is approximately 3 dB. Also, although the unit step 
has the value of zero for positive SJR, the value of 10‘ 5 
has been selected so that the logarithmic scale best shows 
the simulation results (which are all > 10~ 5 ) . As with the 
previous cases examined, there are roughly two regimes of 
operation: interference-dominated and noise-dominated. 
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Figure 58. Probability of Bit Error versus SJR for 
Noncoherent 4FSK with AWGN and Interference. 



Although the curves do not change abruptly from 
decreasing monotonically to remaining constant, approximate 
values that separate the two regions may be identified by 
inspection. For example, the 9-th curve (for SNR=8.99 dB) 
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becomes flat for bit SJR > 12 dB, meaning that the 
performance of the system is noise dominated for bit SJR > 
12 dB. In general, as with previous cases, SNR = SJR is the 
border between the interference-dominated (SNR>SJR) and 
noise-dominated (SJR>SNR) regimes of operation for lower 
values of the error probability. Therefore, for small values 
of the error probability (larger values of SNR and SJR) the 
smaller of the two (SNR or SJR) by and large controls the 
error probability. 

6 . Results For Noncoherent 8FSK 

The probability of bit error for noncoherent 8FSK as a 
function of the signal-to-noise ratio SNR is shown in Figure 
59. Ten curves for signal-to-interference (SJR) ratios from 
-5 dB to +15 dB in 2.223 dB increments are shown. The solid 
line is the probability of error for noncoherent 8FSK 
without co-chanel interference (AWGN only) . We note the 
dramatic increase in the bit error probability due to co- 
channel interference relative to the AWGN only. 
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Figure 59. Probability of Bit Error versus SNR 
Noncoherent 8FSK with AWGN and Interference. 
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The probability of bit error for noncoherent 8FSK as a 
function of the signal-to-interf erence ratio SJR is shown in 



Figure 60. 
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Figure 60. Probability of Bit Error versus SJR 
Noncoherent 8FSK with AWGN and Interference. 
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Ten curves for symbol signal-to-noise (SNR) ratios from 
-5 dB to +15 dB in 2.223 dB increments are shown in Figure 



60. The solid line (step function) shows the probability of 
bit error for the cases when only interference is present 
(no noise) . Since the plots in Figures 58 and 59 represent 
the probability of bit error for noncoherent 8FSK the 'unit 
step' may also be defined on per-bit basis. For M = 8 the 
unit step is shifted to the left by a factor of 

10- log- ( , which is approximately 4.71 dB. Also, although 
log(M) 

the unit step has the value of zero for positive SJR the 
value of 10' 5 has been selected such that a logarithmic 
scale that shows the simulation results (which are all > 
10~ 5 ) can be used. As with the previous cases examined, 
there are roughly two regimes of operation: interference- 

dominated and noise-dominated. Although the curves do not 
change abruptly from decreasing monotonically to being 
horizontal, approximate values that separate the two regions 
may be identified by inspection. For example, the 9-th curve 
(for SNR=8.007 dB) becomes flat for bit SJR > 10 dB meaning 
that the performance of the system is noise dominated for 
bit SJR > 10 dB. In general, as before, SNR=SJR is the 
border between the interference-dominated (SNR>SJR) and 
noise-dominated (SJR>SNR) regimes of operation for lower 
values of the error probability. Therefore, for small values 
of the error probability (larger values of SNR and SJR) the 
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smaller of the two (SNR or SJR) by and large controls the 
bit error probability. 

C . OBSERVATIONS 

The simulation results for the bit error probability as 
a function of SNR with SJR as a parameter can be used to 
determine the increase in bit error probability due to co- 
channel interference (referred here as 'jamming' ) . In the 
case of intentional jamming, the objective is to increase 
the number of bit errors. In such a case, the jammer needs 
to minimize the receiver's signal-to-interference ratio as 
much as possible. The maximum signal-to-interference ratio 
that still allows for the desired increase in the 
probability of bit errors can be determined from the 
simulation results. Given the desired jamming-generated 
increase in bit error probability, the required jammer power 
can be calculated (also given the other relevant parameters 
such as the distance from the target receiver, the receiver 
antenna gain in the direction of the jammer, and the like 
are also known) . For example, from the curves for the 
noncoherent 8FSK case, a -10 dB signal-to-interference ratio 
results in a bit error probability of 0.5, which means that 
the communications is essentially impossible. As SJR 
increases, the bit error probability decreases. 
Nevertheless, bit error probability for both noise and 
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interference is always larger than the bit error probability 
for noise alone. For example, for SJR=10.23 dB and 
SNR=10.23 dB, the bit error probability is about 20 times 
larger than the bit error probability for noise alone. 

In summary, assuming that the operational SNR for the 
MFSK communication channel is known, the family of curves of 
P e versus SNR with SJR as a parameter can be used to 
determine the increase in the probability of error due to 
various levels of co-channel interference. In the case of 
intentional jamming with another synchronized MFSK signal, 
the curves may be used to determine the jammer-on-target 
requirements for a certain increase in bit error 
probability. Subsequently, for a given scenario (jammer-to- 
target distance, etc.) the jammer requirements such as 
jammer power, antenna gain, etc. can be calculated. 

The second set of curves gives the bit error 
probability of MFSK communication systems in the presence of 
both noise and co-channel interference as a function of SJR 
and with SNR ratio as a parameter. The effect of co-channel 
interference is to introduce errors in transmission that are 
quantified by the symbol error probability P e or the Bit 
Error Rate P b (BER) . For BFSK the largest bit error 
probability is 0.5 or 50%. For MFSK the largest symbol error 
probability is given by: 
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Max (P e ) = 1 - 1/M 



(5.1) 



where M is the number of symbols. The maximum bit error 
probability is given by: 



Max (P b )= — ~r) = 0 . 5 (5.2) 

2 M-l M 

Therefore, the maximum bit error probability for MFSK is 
also 0.5, and in the noise-free case (SNR tends to 

infinity) , the theoretical bit error probability due to co- 
channel interference is a unit step function. The curves 
for bit error probabilities when noise and co-channel 
interference are both present tend to 0.5 for negative 
signal-to-interference ratios (in dB) and tend 
asymptotically to zero (no errors) for large positive SJR 
ratios. The curves for negative SNR's are essentially 'flat' 
since the noise 'controls' the errors. On the other hand, if 
the SNR is large and positive (the noise power is small) , 
the curves tend to the 'unit step' function of the 
interference-only case. The most realistic cases of 

comparable SNR and SJR values are those that fall in between 
the 'flat' curves at high (close to 0.5) values of 
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probability of error and the 'unit step' function of the 
interference-only case. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

In this thesis MFSK communication systems with 
coherent and noncoherent detection were simulated in the 
time domain using MATLAB Simulink and the Communications 
Toolbox. The bit error probabilities for coherent and 
noncoherent detection of BFSK, 4FSK, and 8FSK have been 
determined by simulation and verified against the 
theoretical values for the case of AWGN. The MFSK models 
were modified to include the co-channel interference and the 
bit error probabilities obtained for the ranges SNR and SJR 
from -5 dB to +15 dB. 

Simulation results for coherent detection of BFSK, 
4FSK, and 8FSK in the presence of AWGN show excellent 
agreement with theoretical results. On the other hand, 
simulation results for the probability of bit error for 
noncoherent detection of BFSK, 4FSK, and 8FSK differ from 
the theoretical probabilities of bit error by approximately 
-30% for noncoherent BFSK, -20% for noncoherent 4FSK and - 
13% for noncoherent 8FSK, indicating a 'systematic' error in 
the Communications Toolbox implementation of the noncoherent 
MFSK detection. Finally, for symbol signal-to-noise and 
signal-to-interference ratios larger than 15 dB (SNRs and 
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SJR > 15 dB) , the required computational time can be 
excessive for 100 or so errors (days or even weeks on a 200 
MHz PC) . 

B . RECOMMENDATIONS 

The continuation of the research may include: 

• Derivation of the theoretical expression for the 
bit error probability for MFSK with AWGN and co-channel 
interference (BER as a function of SNR and SJR) for the 
following cases: 

• phase-locked co-channel interference 
(relative difference between the interference and signal 
equal 0) , 

• random phase of the co-channel interference 
(relative to the signal) with uniform distribution 
between 0 and 2n, and 

• random phase of the co-channel interference 
with Gaussian distribution. 

• Verification of the theoretical bit error 
probabilities for the above cases using Matlab/Simulink 
models for Monte Carlo type simulations. 
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APPENDIX A. MATLAB PROGRAM SUN_PREP1 .M FOR COHERENT MFSK 
CASE. 

This program prepares the data file for simulation runs 
in case of coherent detection of a MFSK communication system 
in presence of both additive Gaussian noise and co-channel 
interference . 



clear 

noise_only = menu ( ' Select: 

' NOISE' , . . . 

' NOISE and INTERFERENCE ' ) ; 

num_levels = input ('Enter the number of frequencies M [2]: 

' ) ; 

if isempty (num_levels) , num_levels = 2; end 

T_sym = input ('Enter the symbol duration T [1]: '); 

if isempty (T_sym) , T_sym = 1; end 

oversampling = input ('Enter the oversampling factor [2]: '); 

if isempty (oversampling) , oversampling = 2; end 

min_SNR = input ('Enter the MIN Signal to Noise ratio [-5 dB] : 

') ; 

if isempty (min_SNR) , min_SNR = -5; end 

max_SNR = input ('Enter the MAX Signal to Noise ratio [15 
dB] : ' ) ; 

if isempty (max_SNR) , max_SNR = 15; end 
if noise_only ~= 1 

min_SJR = input ('Enter the MIN Signal to Interference 
ratio [-5 dB] : ' ) ; 

if isempty (min_S JR) , min_SJR = -5; end 
max_SJR = input ('Enter the MAX Signal to Interference 
ratio [ 10 dB] : ' ) ; 

if isempty (max_SJR) , max_SJR = 10; end 

end 

if min SNR == max SNR 
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1 ; 



num_noise = 
else 

num_noise = input (' Enter the number of values for SNR 
[ 10 ]: '); 

if isempty (num_noise) , num_noise = 10; end 

end 

if noise_only ~= 1 

if min_SJR == max_SJR 
num_ j am = 1 ; 
else 

num_jam = input ('Enter the number of values for SJR 
[ 10 ]: '); 

if isempty (num_ jam) , num_jam = 10; end 

end 

else 

num_ j am = 1 ; 

end 

min_errors = input ('Enter the min number of errors 
acceptable [100]: '); 

if isempty (min_errors) , min_errors = 100; end 

error_f actor = input ('Enter the factor multiplying the 
number of errors [2]: '); 

if isempty (error_f actor ) , error_factor = 2; end 

initial_num_symbols = error_f actor *min_errors ; 

max_randint = input ('Enter the maximum size of the random 
integer arrays [10 A 6] : '); 

if isempty (max_randint) , max_randint = 10^6; end 

file_name = input ('Enter the file name to save data [no 
ext] : ' , ' s' ) ; 

save sun datal 
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APPENDIX B. MATLAB PROGRAM SUN_FSK1.M FOR COHERENT MFSK CASE 

This program runs coherent MFSK co-channel 

interference with additive Gaussian noise. 



clear 

load sun_datl 
delta_freq = 
f_max = 

delta_t 
seeds = 

signal_seed = 
noise_seed = 
interf seed = 



initial_ 

tic 



%This loads all the input data 

1/T_sym 

= (num_levels-l ) *1/T_sym; 

= 0 . 5/ ( f_max*oversampling) 

= r andint (3,1,1000) ; 
seeds ( 1 ) ; 
seeds (2) ; 
seeds (3) ; 

num symbols = error factor *min errors; 



if num_noise > 1 

delta_SNR = (max_SNR - min_SNR) / (num_noise -1) ; 
else 

delta_SNR = 0; 

end 



SNR = min_SNR + [0:num_noise - 1 ] *delta_SNR; 
noise_var_vect = T_sym/ (2*delta_t ) .* 10 . A (-SNR/10); 

if noise_only ~= 1 
if num_jam > 1 

delta_SJR = (max_SJR - min_SJR) / (num_jam -1) ; 
else 

delta_SJR = 0; 

end 

SJR = min_SJR + [0:num_jam - 1] *delta_SJR; 
interf_gain_vect = 10 . A (-SJR/20) ; 
else 

SJR = - 100; % There is no Jamming so the SJR 

in dB is -infinity 

interf_gain_vect = 0; 

end 



BER = zeros (num_noise, num_j am) ; 
total_symbols = initial_num_symbols; 
for noise_case = l:num_noise 

noise var = noise var vect (noise case); 
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for jam_case = l:nuiti_jam 

enough_errors = 0; 
num_err = 0; 

num_symbols = initial_num_symbols ; 

rand_int = min ( [num_symbols max_randint] ) ; 

total_symbols = num_symbols 

while enough_errors ~= 1 

interf_gain = interf_gain_vect ( jam_case) ; 
clear error_number 

sim ( ' f skco_bm' , num_syiabols) 

[new_errors err_cols] = size (error_number ) ; 
num_err = num_err + new_errors; 
if num_err == 0 

num_symbols = num_symbols*min_errors; 

rand_int = min ( [num_symbols 

max_randint] ) ; 

total_symbols = total_symbols + 

num_s ymb o 1 s ; 

elseif (num_err > 0 & num_err < min_errors) 
num_symbols = (min_errors - 
num_err) *ceil (total_symbols/num_err) ; 

total_symbols = total_symbols + 

num_symbols; 

else 

enough_errors = 1 ; 
end 

end 



number_of_errors (noise_case, jam_case) = num_err 
number_of_symbols (noise_case, jam_case) = 
total_symbols 

BER (noise_case , jam_case) 



= num err / 



total_symbols 












end 












eval ( [ ' 


save ' 


file 


name 


' . snr 


SNR -ascii ' ] ) ; 


eval ( [ ' 


save ' 


file 


name 


' . sjr 


SJR -ascii ' ] ) ; 


eval ( [ ' 


save ' 


file 


name 


' .her 


BER -ascii ' ] ) ; 


eval ( [ ' 


save ' 


file] 


name 


’ .ner 


number of_errors - 


ascii ' ] ) ; 












eval ( [ ' 


save ' 


file_ 


name 


' .nsy 


number of_symbols - 


ascii ' ] ) ; 













end 

toe 
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APPENDIX C. MATLAB PROGRAM SUN_PNC1.M FOR NON-COHERENT MFSK 
CASE 

This program prepares the data file for simulation runs in 
case of non-coherent detection of a MFSK communication 
system in presence of both additive Gaussian noise and co- 
channel interference. 



%%% This prepares the data file for sun runs 

clear 

noise_only = menu ( ' Select: 

' NOISE' , . . . 

' NOISE and INTERFERENCE '); 

num_levels = input ('Enter the number of frequencies M [2]: 

' ) ; 

if isempty (num_levels) , num_levels = 2; end 

T_sym = input ('Enter the symbol duration T [1]: '); 

if isempty (T_sym) , T_sym = 1; end 

oversampling = input ('Enter the oversampling factor [2]: '); 

if isempty (oversampling) , oversampling = 2; end 

min_SNR = input (' Enter the MIN Signal to Noise ratio[-5 dB] : 
' ) ; 

if isempty (min_SNR) , min_SNR = -5; end 

max_SNR = input ('Enter the MAX Signal to Noise ratio [15 
dB] : ' ) ; 

if isempty (max_SNR) , max_SNR = 15; end 
if noise_only ~= 1 

min_SJR = input ('Enter the MIN Signal to Interference 
ratio [-5 dB] : ' ) ; 

if isempty (min_S JR) , min_SJR = -5; end 
max_SJR = input ('Enter the MAX Signal to Interference 
ratio [10 dB] : ' ) ; 

if isempty (max_S JR) , max_SJR = 10; end 

end 

if min SNR == max SNR 
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1 ; 



num_noise = 
else 

num_noise = input ('Enter the number of values for SNR 
[ 10 ]: '); 

if isempty (num_noise) , num_noise = 10; end 

end 

if noise_only ~= 1 

if min_SJR == max_SJR 
num_j am = 1 ; 
else 

num_jam = input ('Enter the number of values for SJR 
[ 10 ]: '); 

if isempty (num_j am) , num_jam = 10; end 

end 

else 

num_ j am = 1 ; 

end 

min_errors = input ('Enter the min number of errors 
acceptable [100]: '); 

if isempty (min_errors) , min_errors = 100; end 

error_f actor = input ('Enter the factor multiplying the 
number of errors [2] : ' ) ; 

if isempty (error_factor) , error_factor = 2; end 

initial_num_symbols = error_f actor*min_errors ; 

max_randint = input ('Enter the maximum size of the random 
integer arrays [ 1 0 A 6 ] : '); 

if isempty (max_randint ) , max_randint = 10^6; end 

file_name = input ('Enter the file name to save data [no 
ext] : ' , ' s ' ) ; 

save sun dncl 
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APPENDIX D . MATLAB PROGRAM SUN_NC1 .M FOR NON-COHERENT MFSK 
CASE 

This program runs coherent MFSK co-channel interference 
with additive Gaussian noise. 

clear 



load sun_dncl %This loads all the input data 
delta_freq = 1/T_sym 

f_max = (num_levels-l) *1/T_sym; 

delta_t = 0 . 5/ (f_max*oversampling) 

seeds = randint (3, 1, 1000) / 

signal_seed = seeds (1); 
noise_seed = seeds (2); 
interf_seed = seeds (3); 

initial num symbols = error factor*min errors; 



tic 



if num_noise > 1 

delta_SNR = (max_SNR - min_SNR) / (num_noise -1) ; 
else 

delta SNR =0; 
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end 



SNR = min_SNR + [O:num_noise - l]*delta SNR; 
noise_var_vect = T_sym/ (2*delta_t) * 10 . A (-SNR/10) ; 

if noise_only ~= 1 
if num_jain > 1 

delta_SJR = (max_SJR - min_SJR) / (num_jam -1); 
else 

delta_SJR = 0; 

end 

S JR = min_SJR + [0:num_jam - 1] *delta_SJR; 
interf_gain_vect = 10 . A (-SJR/20) ; 
else 

SJR = - 100; % There is no Jamming so the SJR 

in dB is -infinity 

interf_gain_vect = 0; 

end 

BER = zeros (num_noise, num_jam) ; 
total_symbols = initial_num_symbols; 
for noise case = lrnum noise 
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noise_var = noise_var_vect (noise_case) ; 

for jam_case = l:num_jam 

enough_errors = 0; 
num_err = 0; 

num_symbols = initial_num_symbols; 

rand_int = min ( [num_symbols max_randint] ) ; 

total_symbols = num_symbols 

while enough_errors ~= 1 

interf_gain = inter f_gain_vect ( jam_case) ; 
clear error_number 
sim ( ' f sknc_bm' , num_symbols) 

[new_errors err_cols] = size (error_n umber ) ; 
num_err = num_err + new_errors; 
if num_err == 0 

num_symbols = num_symbols*min_errors ; 

rand_int = min ( [num_symbols 

max_randint] ) ; 

total_symbols = total_symbols + 

num_symbols ; 

elseif (num_err > 0 & num_err < min_errors) 

num_symbols = (min_errors - 
num_err) *ceil (total_symbols/num_err) ; 
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total_symbols = total_symbols + 

num_symbols ; 

else 

enough_errors = 1; 

end 

end 

number_of_errors (noise_case, jam_case) 

number_of_symbols (noise_case, jam_case ) 
total_symbols 

BER (noise_case, jam_case) 
total_symbols 

end 



eval ( [ ' 


save 1 


1 file name 


1 . snr 


SNR -ascii ' ] ) ; 


eval ( t ' 


save ' 


file name 


1 .sjr 


SJR -ascii ' ] ) ; 


eval ( t ' 


save ' 


file name 


' .ber 


BER -ascii ' ] ) ; 


eval ( [ ' 
ascii ' ] ) ; 


save ' 


file name 


' .ner 


number of errors 


eval ( [ ' 
ascii ' ] ) ; 


save 1 


1 file name 


' .nsy 


number of symbols 



end 



toe 



= num err 



= num err / 
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APPENDIX E. THEORETICAL BER FOR COHERENT BFSK WITH AWGN AND 



CO-CHANNEL INTERFERENCE 

The BER as a function of SNR for negative values of SJR 
(interference power exceeds the signal power) exhibit 
somewhat counter-intuitive behaviour. Namely, the BER 
increases as the SNR increases (there are more errors for 
lower noise power levels) . In order to resolve whether or 
not this is an artifact of the simulation, we have derived 
the expression for the theoretical BER for coherent BFSK 
with AWGN and co-channel interference which is assumed to be 
coherent with the signal. This result of the derivation is: 
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where SJR and SNR are assumed to be in dB. 

Equation (E.l) allows us to plot BER as a surface with 
SNR and SJR as variables, as shown in Figure 60, for the 
ranges of SNR and SJR between -20dB and +20dB. 
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BER as a surface vs SNR and SJR 




P 



error 



Figure 60. Probability of Bit Error versus SNR and SJR for 
Coherent BFSK with AWGN and Interference. 



The probability of bit error as a function of SNR and 
SJR shows a region of shallow local minima with the 
probability of bit error decreasing as the noise power 
increases. This occurs only for high values of interference 
and noise and only up to an SJR of 0 dB. For higher values 
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of SJR there is no local minimum, the decrease with SNR 
becomes monotonic. The conclusion is that noise actually 
reduces the probability of bit error at low SJR. This is 
mainly of academic interest since the reduction occurs at 
very high values of error probability (from about 0.5 to 
about 0.35) . 

Equation (E.l) can also be used to generate conventional 
BER plots. These plots represent cross-sections of the 
surface illustrated in Figure 60 along planes of constant 
SNR, shown in Figure 61, or constant SJR, shown in Figure 
62. 




,.-20, SNR ,20, 

*■ m 

Figure 61. Probability of Bit Error versus SNR for 
Coherent BFSK with AWGN and Interference. 
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Six curves for SJR ratios from -20 dB to 0.4 dB in 4.1 



dB increments are shown in Figure 61. 

The theoretical BER for coherent BFSK as a function of 
the SJR is shown in Figure 62. Seven curves for SNR ratios 
from -20 dB to 4 . 5 dB in 4.1 dB increments are shown. 




u -20j SJR .20.J 

n 

Noise Variance 

Figure 62. Probability of Bit Error versus SJR for 
Coherent BFSK with AWGN. 
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